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ABSTRACT 



The deconfinement transition in SU(4) lattice gauge theory is investigated 
on N% x N t lattices for N s = 8-16 and N t = 4-8 using a modified Wilson 
action which is expected to be free of any bulk transitions. The susceptibility 
X|l| X ' where L is the order parameter for deconfinement, is found to increase 
linearly with spatial volume for N t = 4, 5, and 6, indicating a first order 
deconfinement phase transition. The latent heat of the transition is estimated 
to be ~ | of the corresponding ideal gas energy density at T c . 
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1. INTRODUCTION 



It is widely expected that the ongoing experiments at the Relativistic 
Heavy Ion Collider, BNL, New York, and those at the proposed Large Hadron 
Collider, CERN, Geneva, will provide us a strong evidence for Quark-Gluon 
Plasma and may even yield crucial details about this new phase and the 
nature of the phase transition(s) leading to it. The physics driving this 
phase transition is naturally of great interest. While the real world has 
two light (u,d) flavours of quarks and one somewhat heavier (s) flavour, 
both analytical and numerical methods in lattice Quantum Chromodynamics 
(QCD) begin from the limiting cases of either massless or infinitely massive 
quarks. One talks of the chiral symmetry restoring phase transition and 
the deconfinement phase transition in these two cases respectively and has 
suitable order parameters to investigate them. For quarks with N colours 
and Nf flavours, these transitions are related to spontaneous breaking of a 
global Z(N) and SU(Nf) x SU(Nf) chiral symmetry, which QCD has in the 
infinite quark mass and zero quark mass limit respectively. Which of them is 
more relevant in the real world, is a priori not clear, since these symmetries 
are broken explicitly to various extents. The low masses of the light flavours 
suggest chiral symmetry to be the dominant one, leading one to expect a 
chiral transition in QCD at finite temperature. Indeed, the behaviour of the 
corresponding order parameter, the chiral condensate (ipip), confirms this 
expectation in numerical simulations. However, it is also seen in the same 
numerical simulations that even the order parameter for the deconfinement 
transition, 



acquires nonzero values at this chiral transition; in principle, it could have 
done so already at temperatures a lot below the chiral transition since no 
symmetry prevents it from doing so. The energy density also shows a large 
change at the transition which would be unexpected from a naive count 
of degrees of freedom at a chiral transition but is in fact consistent with 
that expected of a deconfinement transition. These apparently mysterious 
observations for QCD with light dynamical quarks can be explained |I] using 
arguments based on the large N limit, if the deconfinement transition for 
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N > 4 is of second order. SU(4) is clearly the simplest case to check whether 
this is so. 
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Figure 1: A schematic phase diagram in (/3, Pa) -plane for the mixed action 
ofeq. (g). 



Numerical simulations of SU(4) theory at finite temperatures have been 
done in the past p| and recently || as well. All of them reported a first order 
deconfinement transition but they all used the Wilson action, or the more 
general mixed action: 

Here Up denotes the directed product around an elementary plaquette P 
of the link (gauge) variables, U^x), located at site x in the direction \i 
=1-4, tr and trA denote the traces of Up in the fundamental and adjoint 
representations, the /3's are the corresponding couplings, and N is the number 
of colours. The sum over P is over all independent plaquettes. A well known 
problem in the simulations with these actions, especially for large N, is the 
presence of bulk transitions which are lattice artifacts. The generic phase 
diagram of the mixed action (||D in its coupling plane is depicted schematically 
in Fig. The solid lines in it show first order bulk transition lines. The 
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dotted line after the end point D is drawn to suggest the impact D may 
have on the Wilson axis {(3a = 0). For iV > 4, D is expected to be where 
E is shown, causing a first order bulk transition for the usual Wilson action. 
In order to avoid it, simulations were made at negative (3a @ and/or for 
larger 1^ = 6 and obtained a first order deconfinement phase transition 
for SU(4). 

Extensive studies [|J of the deconfinement phase transition for the action 
above for the N = 2 case have, however, shown that the presence of bulk 
transitions affect the order and location of the deconfinement transition in 
subtle and somewhat inexplicable ways. For example, recall that one expects 
the coupling at which the deconfinement transition takes place to move log- 
arithmically with N t , the temporal extension of the lattice whereas a shift 
in the bulk transition point is expected to be inversely proportional to the 
4-volume of the lattice, N^N t , where N s is the spatial lattice size. This was 
the rationale behind the choice of simulation parameters in the earlier [^]. [| 
studies of the 577(4) theory. However, we [f§] found that either this expec- 
tation is incorrect or there are apparent qualitative violations of universality 
for the SU (2) theory since its deconfinement transition was of second order 
for small (3a but was of first order for large enough (3a for an entire range of 
N t from 2 up to 8, with hardly any shift with N t . Indeed, universality was 
restored || in that case only after eliminating the bulk transitions associated 
with the Z{2) vortices and Z{2) monopoles by adding large chemical poten- 
tials for them. It seems natural to expect that the bulk transitions for N = 4 
may also have affected the earlier results on the order of its deconfinement 
transition for the action above for Nt — 4 or 6. A cleaner determination of 
the order of the deconfinement transition for S77(4) may be obtained by wip- 
ing off the corresponding bulk transitions by suppressing the relevant Z{N) 
objects. We pursue this idea here to investigate the deconfinement phase 
transition in the 577(4) theory. 

2. SIMULATIONS AND RESULTS 

Generalizing the idea of positive plaquette models || in the literature for 
the SU{2) lattice gauge theory to the case of SU (N) theories, we propose to 
simulate the action 

^(i-^M-M), (3) 
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for N = 4, where —it < a < it is the phase of tr Up. By adding the adjoint 
term of eq. (2) to the action (3), one sees that the phase diagram of the 
resultant mixed action should not have the bulk lines AC or BC and hence 
the endpoint D or E. This is because the action (|]) restricts tr Up to lie only 
in one Z(iV)-sector for all, even small, (3. 

We have simulated the above action on x N t lattices for N s =8, 10, 
12, 15, 16 and N t = 4, 5, 6, 8 using a 15-hit Metropolis et al. algorithm. 
The calculations were done on a cluster of pentiums. Typically short runs to 
look for points of rapid variations in (\L\) were followed by long runs (a few 
million sweeps) to determine the susceptibility x\l\, 

X\l\ = N! «|L| 2 ) - <|L|> 2 ) , (4) 

as a function of f3 using the histogramming technique @ . In order to monitor 
the presence of any bulk transition, we also studied the plaquette suscepti- 
bility given by, 

Xp = GN*N t ((P 2 > - <P> 2 ) , (5) 

where P denotes the average of (Re tr £7p)/N over the entire lattice. Usual fi- 
nite size scaling techniques were used to determine the order of the transition 
and its exponents. According to the finite size scaling theory ||, the peak of 
the \L\ (or plaquette) susceptibility at the location of the deconfinement (or 
bulk) transition should grow on N% x N t lattices like 

X^for P OC K , (6) 

for fixed Nt. For a second order transition, u = j/u, where 7 and v char- 
acterize the growth of the \L\ (plaquette)-susceptibility and the correlation 
length near the critical temperature (coupling) on an infinite spatial lattice. 
If the phase transition were to be of first order instead, then one expects 
the exponent u = 3, corresponding to the dimensionality of the space. In 
addition, of course, the (\L\) or (P) is expected to exhibit a sharp, or even 
discontinuous, jump. The corresponding probability distribution should then 
show a double peak structure in stead of a usual gaussian distribution. 

In simulations on iV s 3 x 4 lattices, N s = 8, 10, 12, we found the hot (all 
link variables random) and cold (all links equal to identity) starts to converge 
quickly to a unique value of \L\ at couplings a little away from the transition 
point on its both sides but a clear co-existence of states was visible for all 
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Figure 2: The evolution of \L\ as a function of Monte Carlo time for x 4 
lattices close to the transition point. 
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Figure 3: The histograms of \L\ on iV s 3 x 4 lattices for long runs close to the 
transition point. 
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Figure 4: The susceptibility as a function of /3 for iV 3 x 4 lattices. 
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Figure 5: The histograms of the plaquette P on N% x 4 lattices for the same 
runs as in Fig. |[ 

lattices at the transition point, as shown in Fig. |2]. It is clearly seen that 
the tunnelling frequency goes down with spatial volume. The histograms of 
\L\, shown in Fig. |3], display peaks which become narrower with increasing 
volume while the gap between them remains unchanged. These classic signs 
of a first order phase transition are confirmed by a quantitative analysis of 
the growth of with volume, as seen in Fig. |j. The horizontal lines in 
each case are predictions obtained by scaling the N s = 8 results linearly with 
volume, as expected from eq. (Q) and u = 3 for a first order transition. 

Although the above first order transition is clearly a deconfinement tran- 
sition, as the behaviour of the order parameter \L\ certifies, it could, of course 
be due to a coincident bulk transition. One can check for this possibility by 
analogous studies of the average plaquette P. Fig. | displays the histograms 
of P for the same runs as in Figs. |2] and [| These histograms do exhibit 
curious non-gaussian structures, with a suggestion to develop a discontinuity 
as the spatial volume grows. The corresponding susceptibilities turn out to 
be too noisy to be conclusive, especially when compared to the behaviour of 
|L|-susceptibility in Fig. [|. However, their growth with volume is not incon- 
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sistent with being linear. We therefore made simulations on larger symmetric 
lattices up to the sizes 12 4 at these couplings and in their immediate neigh- 
bourhood and found only usual gaussian distributions for the plaquette. For 
a genuine bulk transition, the behaviour in Fig. |5| should have been accentu- 
ated with the three-fold increase in the 4-volume. Recall that the presence 
of latent heat at a first order deconfinement phase transition can also be a 
source of the behaviour in Fig. [|. This should then lead to specific predic- 
tions as the temporal lattice size is increased. We postpone discussing them 
in the next section, turning in stead to investigations of larger N t . 

For larger Nt, we used many longer runs in the region of strong variation 
of (\L\) to obtain the susceptibility directly and used the histogramming 
technique only for the finer determinations of the critical coupling. Our 
results for (|L|) as a function of (3 clearly show the expected shift for a 
deconfinement phase transition for N t = 4 and 6, as can be seen in Fig. || 
As is usual for the SU (2) and SU (3) theories, the rise of the order parameter 
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Figure 6: \L\ as a function of f3 for lattices with N t 4 and 6. 



at p c is seen to be slower in Fig. || as N t is increased due to the well-known 
ultra-violet effects. The j3 c values, determined from the xjzf* an d listed in 
Table [I] for N t = 4, 5 and 6 for different spatial volumes, also evidently 
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Table 1: The values of (3 at which long simulations were performed on iV 3 x N t 
lattices, (3 C and the height of the |L (-susceptibility peak, xfzT- 



N* x N t 


P 


(3c,N s 


.max 
A. \L\ 


v max 

A predicted 


8 3 x 4 
10 3 x 4 
12 3 x 4 


10.360 
10.364 
10.363 


10.360(2) 
10.364(5) 
10.363(5) 


5.48(16) 

11.16(35) 

22.41(1.11) 


10.71(32) 
18.5(6) 


10 3 x 5 
15 3 x 5 


10.520 
10.520 


10.515(5) 
10.525(5) 


4.63(40) 
14.02(1.50) 


15.6(1.3) 


12 3 x 6 
16 3 x 6 


10.675 
10.675 


10.686(5) 
10.676(5) 


4.36(35) 
10.43(95) 


10.3(8) 



suggest strong shifts in j3 c with N t . Comparing our results in Table |T] with 
those of Ref. for N t = 6 for the usual Wilson action, i.e, for eq.(H) with (3a 
= 0, one finds a shift by ~ 0.1 towards smaller coupling. Using in each case 
the peak height for the smaller spatial volume, the ^ max on the bigger lattice 
can be predicted, assuming a first order deconfinement phase transition. The 
predictions listed in Table |^ can be seen to be in very good agreement with the 
direct Monte Carlo determinations. Along with the shifts in /3, these confirm 
that the same physical first order deconfinement phase transition is being 
simulated on these lattices as the continuum limit of a — > is approached in 
a progressive manner by increasing the temporal lattice size N t . 



3. SCALING AND LATENT HEAT 



A quantitative test of the fact that the first order transition reported 
in the previous section is indeed a physical deconfinement transition (and 
not a bulk transition) consists of translating all the (3 C in Table [j] to the 
corresponding transition temperature in physical units and thus checking 
whether it is constant while the /3 c 's shift with Nt as expected from the 
renormalization group equation, 

_, 1 T c /46„r'""» / B \ ... 

= a! \ p) exp bd ■ (7) 
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Figure 7: 7V t 1 as a function of /3 along with the RG-curve mentioned in the 
text. 
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Figure 8: Nf 1 as a function of 0e along with the corresponding RG-curve. 
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where 

60 = 12^' and 61 = ' (8) 

are the first two coefficients of the perturbative /3-function for the S77(4) 
Yang-Mills theory. Fig. [7| exhibits a comparison of our set of /3 C for N t =4, 
5, 6 and 8 with eq.(^). The normalization was chosen by demanding the line 
to pass through the point at Nt = 6. While the shifts are in accord with 
the expectations, one also sees strong quantitative deviations. These are, 
however, not unusual: similar deviations have been seen in the studies for N 
2 and 3 |T(J as well. As pointed earlier [IT|, and tested successfully [|10j| 



for SU{2) and SU(3), one possible cure for removing these deviations is to 
employ a better choice of the coupling in eq.(^) which consists of replacing 
(3 there by 

_ 0-25(iV 2 - 1) 
fa ~ 1 _ {p) ■ (9) 

Fig. |8| displays our data using this new variable along with the corresponding 
RG-curve of eq.(^) normalized the same way as before. One sees an excel- 
lent agreement, providing a concrete quantitative evidence in favour of the 
physical nature of the transition. Furthermore, it leads to an estimate of T c 
= 18.5 ± 0.5Ai, where the error is estimated by requiring all the data points 
to lie within the band generated by it. 

Apart from the characteristic (logarithmic) shift of the transition point 
with N t , seen above, the latent heat of a first order deconfinement phase 
transition provides yet another quantitative test of the continuum limit since 
it should also remain constant as Nt — * 00. Requiring the pressure to be 
continuous at the deconfinement phase transition, the latent heat can be 
obtained from two different observables Ai = A(e — 3p)/T c 4 , and A2 = 
A(e +p)/T c 4 , where A denotes discontinuities across the transition in the 
respective variables, 



Ai = -48A, 4 a^— AP 



dg- 2 
> U -da 



A 2 = 32A t 4 ^p(AP t - AP S ) , (10) 

and C(g 2 ) = (1 - 0.2366# 2 + 0{g 4 )) for SU(4) §. We will employ the 
perturbative /5-function in Ai with its coefficients given by eq. (§). In order 



12 




13 



Table 2: Both the latent heat estimates of eq.(|T0|) as a function of N t . 



N t 


4 


5 


6 8 


Ai 


21.03(5) 


11.02(6) 


8.31(5) 6.57(16) 


A 2 


9.89(14) 


7.77(40) 


6.04(60) 6.45(99) 



to obtain the AP, AP S and AP t , the minimum of the histogram iV(|L|) was 
used to separate the two phases in each case. The errors were estimated 
by varying it within the limits suggested by the corresponding histograms. 
From eq. fllOD , it is clear that the plaquette discontinuity AP oc Nf 4 if one is 
to obtain the same latent heat in physical units, since the rest of the factors 
vary mildly with Nt- In case of a bulk phase transition, on the other hand, 
one would expect the plaquette discontinuity to remain constant or exhibit 
a mild increase. Indeed, as Fig. |9] displays, its decrease with N t is seen to be 
consistent with these expectations for N t > 5. Our estimates of latent heat 
from Ai and A 2 are shown in Fig. [10] as a function of N t (filled symbols) 



along with the corresponding results of Ref. fU for N t = 6 (open symbols) 
which can be seen to be in good agreement with ours. Both the estimates in 



eq.pOf) must agree with each other as Nt — > oo, since the neglected cut-off 
corrections become then insignificant. Table and Fig. |H] verify this to 
be the case. One may therefore take the results for N t =8 as a reasonable 
approximation to the continuum limit. Comparing them with the energy 
density of the corresponding ideal gas at T c , €sb = 7r 2 (iV 2 — 1)T 6 4 /15 = vr 2 T c 4 , 
one obtains for the latent heat, 



— = 0.6657 ±0.0162 . (11) 
Here we have used Ai to get the result above. Using A 2 leads to 0.65 ± 0.10. 
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4. SUMMARY AND DISCUSSION 



Establishing the order of the deconfinement phase transition for the £77(4) 
lattice gauge theory at finite temperature is important for our understanding 
of the physics of the phase transition to quark-gluon plasma. Possible subtle 
influences, which the bulk transitions for the mixed action of eq. @ (see 
Fig. |l|) may have, necessitate an approach to bypass them. In this paper, we 
have attempted this by generalizing the idea of positive plaquette model for 
the £77(2) case and simulating the action of eq. (|3|) to investigate the £77(4) 
theory at T ^ on iV s 3 x N t lattice with iV s > 2N t and N t = 4, 5, 6 and 8. 

Various qualitative indicators, such as, the histograms and evolutions of 
the order parameter L, suggest a first order deconfinement phase transition 
for SU(4) on N t = 4 lattices with N s varying up to 12. The linear growth of 
X'\l\ k w ith volume for N t = 4, quantitatively confirms this finding. Increasing 
N t to 5 and 6, one again finds a growth in xfjT' that is consistent with being 
linear in volume. Defining the transition coupling as the location at which 
the maximum of the above susceptibility occurs, one finds a significant shift 
in it as N t is varied from 4 to 8. The amount of shift is consistent with the 
expectations of a physical deconfinement transition; a bulk transition would 
have had much smaller shifts. Indeed, the set of (3 C for N t =4, 5, 6 and 
8 is consistent with the asymptotic scaling relation of eq.(^) when used in 
conjunction with the improved coupling of eq. (|]), leading to a determination 
of the transition temperature T c = (18.5 ± 0.5)Al. 

The presence of latent heat at a first order phase transition necessarily 
means a discontinuity in the average plaquette on a finite lattice, which can 
however also be due to an unphysical bulk transition. We find that the 
plaquette discontinuity AP decreases with increase in N t as approximately 
the fourth power of N t . This indicates both a lack of a bulk transition and 
the presence of a first order deconfinement phase transition, since it suggests 
i) a vanishing discontinuity in the N t —>■ oo limit ii) but in a manner that 
leaves the latent heat constant in physical units. We estimate the latent 
heat to be (0.666 ± 0.016) times the corresponding ideal gas energy density 
at T c using the perturbative /5-function in eq. fllPl). This value, obtained 
from our Nt =8 simulation, is expected to be close to the continuum limit 
since the two different estimates from eq. ( |i~0|) coincide. Comparing with 
the corresponding value (using again perturbative /^-function) for the SU (3) 
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theory, which is about 0.43 [E|, one finds that the deconfinement transition 
grows stronger in nature as the number of colours N is increased. Whether 
this continues as N is increased further is at present not clear. However, 
already the known results for iV = 3 and 4 are sufficient to cast doubt on the 
explanation proposed in Ref . |] . Consequently, the mystery of the behaviour 
of the energy density and the deconfinement order parameter, L, mentioned 
in the introduction remains. 



5. ACKNOWLEDGMENTS 



It is a pleasure to acknowledge interesting discussions with Sourendu Gupta. 



16 



References 



[1] R. D. Pisarski and M. Tytgat, |hep-ph/9702340| . 

[2] G. G. Batrouni and B. Svetitsky, Phys. Rev. Lett. 52 (1984) 2205; 
A. Gocksch and M. Okawa, ibid, 52 (1984) 1751. 

[3] M. Wingate and S. Ohta, Phys. Rev. D63 (2001) 094502. 

[4] R. V. Gavai and M. Mathur, Phys. Rev. D 56 (1997) 32; M. 
Mathur, R. V. Gavai, Nucl. Phys. B448 (1995) 399; R. V. Gavai, 
M. Grady and M. Mathur, Nucl. Phys. B423 (1994) 123. 

[5] R. V. Gavai and M. Mathur, Phys. Lett. B458 (1999) 331; S. 
Datta and R. V. Gavai, Phys. Rev. D 62 (2000) 054512. 

[6] G. Mack and E. Pietarinen, Nucl. Phys. B205 [FS5] (1982) 141. 

[7] A. M. Ferrenberg and R. H. Swendsen, Phys. Rev. Lett. 61 (1988) 
2635. 

[8] M. N. Barber, in Phase Transitions and Critical Phenomena, vol. 
8, Ed. C. Domb and J. L. Lebowitz (Academic Press, New York, 
1983) p. 146. 

[9] M. S. Challa, D. P. Landau and K. Binder, Phys. Rev. B34 
(1986) 1841. 

[10] J. Fingberg, U. Heller and F. Karsch, Nucl. Phys. B392 (1993) 
493. 

[11] G. Parisi, Proceedings of the XXth Conference on High Energy 
Physics, Madison, 1980. 

[12] B. Beinlich, F. Karsch and A. Peikert, Phys. Lett B390 (1997) 
268. 



17 



